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Abstract 



We present experimental evidence and theoretical arguments showing that the time- 
evolution of freely decaying 2-d turbulence is governed by a discrete time scale invariance 
rather than a continuous time scale invariance. Physically, this reflects that the time- 
i-rt ■ evolution of the merging of vortices is not smooth but punctuated, leading to a prefered 

P^ ' scale factor and as a consequence to log-periodic oscillations. From a thorough analysis of 

O ■ freely decaying 2-d turbulence experiments, we show that the number of vortices, their ra- 

dius and separation display log-periodic oscillations as a function of time with an average 
log-frequency of « 4 — 5 corresponding to a prefered scaling ratio of ~ 1.2 — 1.3. 



1 Introduction 

Hydrodynamical turbulence is one of the major remaining challenges in the physical sciences. 
It provides a fascinating phenomenology with both simple universal scaling laws and complex 
system specific features. Two-dimensional turbulence is of special interest, both for its applica- 
tions in astrophysics and geophysics and its theoretical properties. It is also more amenable to 
detailed experimental, numerical and theoretical studies. Its main characteristic property is the 
formation of and interaction between coherent structures or vortices. 

Kolmogorov first proposed his energy cascade model which was later refined by many authors 
to account for intermittency effects (see for instance [|l|). It remains to be assessed to what 
extend this class of statistical phenomenological models is physically correct and whether the 
the cascade picture can be taken to represent real physical processes that occur genuinely in 
fluids. 

2-d turbulence is distinct from 3-d turbulence by the conservation of vorticity along fluid 
particle paths, when ignoring viscosity and forcing (there are no vortex-stretching term in the 
2-d equations). An inverse energy cascade was conjectured by Kraichnan [Q] and has since been 
largely confirmed by simulations and experiments (see for instance Ref . |||] and references therein) . 
This inverse cascade does not appear to be characterised by the presence of strong coherent 
vortices. Similarly, if the direct cascade exists in 3-d turbulence, vortex filaments are also rare 
events but, according to the so-called hierarchical structure model Q, they nevertheless play 
an important role. These coherent structures are the most singular, the most intermittent, and 
the rarest events in a turbulent medium; these structures are also the most sensitive structures 
depending on the environment. According to this theory, "softer" fluctuations are organised 
statistically into a hierarchy which satisfies universal statistical symmetry with scaling exponents 
controlled by the nature of these most singular structures. 

Returning to the 2-d inverse energy cascade, the coherent vortices are clearly identified at 
the forcing scale, comparable to the system size. Their statistical properties suggest that the 
cascade is driven by a clustering mechanism involving vortices of the same sign Q , in contrast to 
the sequence of merging events producing larger and larger vortices occurring in freely decaying 
2-d turbulence which we analyse below. In the inverse cascade, the energy is thus transported 
to larger scales by an aggregation of vortices of the same sign while no apparent fusion is 
present, in contrast to freely decaying 2-d turbulence. In addition, the inverse cascade appears 
non-intermittent with regular Gaussian statistical properties of the velocity increments. The 
generalisation of our finding obtained in a situation of vortex merging to forced 2-d turbulence 
where vortices do not merge but aggregate and its relevance to the energy cascade (inverse in 
2-d and direct in 3-d) are thus conjectural and remain to be tested directly. 

Here, we identify a novel signature of 2-d freely decaying turbulent flows suggesting that 
vortex merging occur via a cascade process with a discrete hierarchy of sizes. We document log- 
periodic oscillations, which are the hallmarks of a discrete hierarchical structure with a prefered 
scaling ratio, in the time evolution of the number of vortices, their radius and separation. 

The theory of log-periodicity (and its associated complex exponents) has advanced signifi- 
cantly Ip in the last few years. Log-periodicity reflects a discrete scale invariance, i.e., the fact 
that dilational symmetry occurs only under magnification under special factors, which are arbi- 
trary powers A" of a prefered scaling ratio A. Log-periodicity has been studied in the eighties in 
relation to various problems of physics embedded in discrete hierarchical systems. In the context 
of turbulence, shell models construct explicitly a discrete scale invariant set of equations whose 
solutions are marred by unwanted log-periodicities. Only recently has it been realized that 
discrete scale invariance and its associated complex exponents may appear "spontaneously" in 
Euclidean systems, i.e., without the need for a pre-existing hierarchy (see P] for references and 



discussion of Laplacian growth models, rupture in heterogeneous systems, earthquakes, animals, 
financial crashes among many other systems). In addition, general field theoretical arguments 
[y] indicate that complex exponents are to be expected generically for out-of-equilibrium and 
quenched disordered systems. 

The paper is organised as follows. In section 2, we discuss the plausibility of log-periodicity 
in turbulence from the vantage of three theoretical arguments based on similarity analysis, 
mean field theory and the aggregation analogy. In section 3 we describe our statistical tools, 
in particular the Lomb periodogram technique and its performance on synthetic tests with a 
detailed discussion of the role of noise and the problem of averaging. These considerations 
concerning detection of log-periodic signals are offered to point-out important pit-falls, which 
might explain why such signals have not previously been reported. The experimental set-up, 
main experimental results and our analysis of the experimental data are presented in section 4 
while section 5 concludes. 

2 Theory 

There is no established theory on 2-d freely decaying turbulence. All available approaches 
exhibit some weaknesses or incompatibilities with experiments. Here we discuss three of the 
main theoretical approaches only to suggest that log-periodicity is not unreasonable from a 
theoretical view point. 

2.1 Dimensional analysis and Similarity of the third kind 

Scaling laws are tools which can be used to quantify the complexity of turbulent flows. Well- 
known examples include the Kolmogorov energy-cascade from large scales to smaller scales in 
3-d turbulence and from small scales to larger scales in forced 2-d turbulence [Q], as already 
mentioned. 

For the decaying 2-d case, the flow carries coherent vortices Q. Let us assume that the 
number of vortices also obeys a scaling law. Using the 11 theorem of dimensional analysis [^, 
which states that the laws of physics must be such as to be expressed using only dimensionless 
variables, it can be written as 

n{t)^^F{Re,^,^,...) . (1) 

Here E is the kinetic energy per area, Re is the Reynolds number inversely proportional to the 
viscosity v, t^ = L/{v) is the typical time scale for a vortex to move over the system size L and 
r = {R)/{v) is the typical time for a vortex to move over a distance comparable to its typical 
radius. That the first term of the r.h.s. is inversely proportional to Et"^ stems from the fact that 
E is proportional to the inverse of a time-squared. 

The simplest assumption is that the function F goes to a constant when the viscosity 
z^ — > 0, leading to the prediction n{t) ~ p-. This corresponds to the so-called complete similarity 
of the first kind S with respect to the variables Re, t/ti and t/r. 

The existence of a finite limit of F was first questioned by L.D. Landau and A.M. Obukhov 
in the context of fully developed turbulence in 3-d, on the basis of the existence of intermittency 
- large fluctuations of the energy dissipation rate about its mean value. Following the same 
argument and Barenblatt's classification leads to the possibility of an incomplete similarity in 
the variables t/ti and t/r. This would require the absence of a finite limit for F{Re, -^, ^, ...) 
as t/tL ^ or as t/r -^ oo, and leads in the simplest case to the form 

n(i) ~ ;^ (^r (^)^ ~ ^ , withe = 2-a-/3, (2) 



where a and /? are intermittency exponents. If a and /3 are real, this corresponds to a similarity 
of the second kind [Q . 

Incomplete self-similarity ||l^, ^ may stem from a possible dependence of the exponents on 
Re. Dubrulle [12| (see also ||l^) has proposed to call the case where a and/or /? is complex, 



leading to 

n(t)~^ cos[e/lnt + V] , (3) 

a similarity of the third kind, characterised by the absence of limit for F{Re, -^, ^, ...) and 
accelerated (log-periodic) oscillations. 

To our knowledge, Novikov has been the first to point out in 1966 that structure functions 
in 3-d turbulence should contain log-periodic oscillations [Q|. His argument was that if an 
unstable eddy in a turbulent flow typically breaks up into two or three smaller eddies, but not 
into 10 or 20 eddies, then one can suspect the existence of a prefered scale factor, hence the 
log-periodic oscillations. This view does not seem to be born out by direct observation of 3-d 
flows where vortex filaments are very rare events and the cascade seems to occur via a different 
mechanism. The dynamical occurrence of a discrete hierarchy of velocity fluctuation remains 
to be investigated. Notice however that our finding of a discrete hierarchy of merging events in 
freely decaying 2-d turbulence is in a sense reverse to Novikov's picture. 

Log-periodic oscillations have been repeatedly observed in structure factors but do not seem 
to be stable and depend on the nature of the global geometry of the flow and recirculation 
|||, |l^ as well as the analysing procedure. As to be discussed in section ^, the detection of 
log-periodic structures in data is a non-trivial task and is rather sensitive to the implementation 
of the analysis. 

2.2 Mean field theory 

At times sufficiently large so that the individual mean time between two collisions suffered by 
any one vortex is larger than the time needed for a typical vortex to move the distance of its 
typical diameter, we can use Trizac's mean field approach [^] based on dimensional analysis. 
The fundamental equation resulting from this analysis governs the time evolution of the vortex 
density n{t) : 

^ = -e^ (4) 

dt ^ t ' ^ ^ 

It integration gives the usual power law solution 

n{t) ~ t-« , (5) 

where ^ is empirically found equal to 0.72 it 0.03 while theoretical arguments suggest ^ = 1 



|13, y] or ^ = 0.8 |6| 



It is important to recognise that (^) represents an averaging over highly intermittent and 
"non-smooth" processes : most of the time, the dynamics of the freely decaying 2-d turbulence 
is dominated by the mutual advection of the vortices, while rare dissipative events of energetic 
vortex merging punctuate this dynamics. Thus, in reality, the time dynamics is not scale free 
as expression (^) would lead us to believe. There are rather well-defined instants at which the 
coalescence proceeds, corresponding to a discrete time evolution for the number of vortices, and 
the occurrence of theses instants have to obey a scale invariance in time dilation in order to 
recover (Q) on average. 

Let us now go beyond this average description and attempt to describe this punctuated 
dynamics by introducing time scales in the simplest way such as to recover the continuous time 
scale invariance (CTSI) upon averaging : we hypothesise that the sudden coalescence can be 



accounted for by breaking partially the CTSI captured by expression (^) into a weaker discrete 
time scale invariance (DTSI). Quantitatively, this amounts to replacing (Q) by expressing dt as 
a discrete time step approximation with should also be scale invariant : 

dt = at , (6) 

with Q a characteristic scale ratio. Then, (Q) transforms into 

n(t + at) — n{t) ^n{t) 

7 = ~^~r~ ' 

at t 

leading to 

n{Xt) = (1 - (A - 1)0 n{t) , where A = 1 + a . (7) 

Expression (0) is now a discrete scale invariant relation for n{t). Its general solutions are power 

laws t~*™ like (^), where the exponents Sm now depend on a and can acquire an imaginary part. 

It is this imaginary part which is the signature of a discrete hierarchy of time scales reflected in 

log-periodic oscillations : 

ln(l - (A - 1)£) 2tt 

'■" = - inA +""hrA- <«> 

where i = \f—\ and m is an integer. The real solution sq recovers the value ^ in the limit 
a — > 0, A ^ 1, where the discrete scale symmetry recovers the continuous scale symmetry. 

We see that, for A > 1, sq > i^. As an example, take the experimentally determined value 
i = 0.72, which gives sq = 0.80 for A = 1.13 and sq = 1-01 for A = 1.4. For A = 1.2, we get 
So = 0.85. This suggests an explanation for the difference between prediction from theory and 
experimental and numerical data : the measurements that assume CTSI underestimate the true 
value that is impregnated by the existence the punctuated dynamics leading to DTSI. 

In this formalism, one can take into account the presence of fluctuations in the rate of vortex 
mergers by replacing (0) by 

n(Aito) = (1 - (Ai - 1)6) n(to) (9) 

n(A2Aito) = (1 - (A2 - 1)6) ^(Aito) (10) 

(11) 

N N 

n(J[ Xito) = n{to) n(l - (^i - 1)^0 • (12) 

1=1 i=l 

These equations describe the fact that the hierarchy of time scales at which the vortex mergers 
occur do not need to obey an exact hierarchy but can exhibit fluctuations. Looking for a power 
law behaviour n{t) ~ t~^, we find that the typical value for s is given by Q 

(ln(l-(A-l)6) , . 27r 

'- = (MX) + ^"^(h^' ('') 

where the averaging (.) is performed over the distributions of A and 6 

2.3 Linear theory of aggregation 

Benzi et al. [f9l and McWilliams jEH] have noted the analogy between the vortex merging 



in turbulence and the kinetic theory of colloidal aggregation |21[| . Starting from a population 



ni{t = 0) of monomers, the kinetic equations of aggregation read 

oo 
-711 ^ Ji,j Uj (14) 



dni 

~dr 



drij 
IF 



i=l 
j-1 oo 



The first equation expresses that the population of monomers decreases by coalescence with all 
possible species. The second equation expresses the competition between the formation of a 
j'-mer by coalescence of an i-mer and a j — z-mer and its destruction by coalescence with all 
possible species. The Jij express the rates at which these coagulation occur. When diffusion 
motion dominates, we have 

Jijr^{Ri + Rj){Di + Dj) , (16) 

where Ri (resp. Di) is the radius (resp. diffusion coefficient) of the i-mer. From Stokes law and 
Einstein's relationship, D ~ 1/R the kernel are homogeneous: J\i,\j = Jij- In this case, the 



-Anin (17) 



Jij can be taken constant and ( p^ ) reduces to 

dni 
If ~ 

dn A 2 '^ 2 MQA 

— - = — ni n , (18) 

where n{t) = Yl^i '^j is the total population. It is then straightforward to check that ni{t) ~ 
n{t) ~ 1/t at large t. 

A fundamental physical ingredient in this approach is the quadratic pair-wise nature of the 



interactions described by the r.h.s. of the equations (15,18), which embodies the idea that 
coagulation is a stochastic process requiring the chancy collision of two particles that can come 
from arbitrary position in space. We argue quite reasonably that this is not the correct situation 
for vortex fusion in freely decaying 2-d turbulence due to the confinement that a 2-d space 
represents. Experiments (see for example figure 3 of Ref . ||22[| ) show that, between two merging 
events, vortices move as in a flexible "cage" of neighbouring vortices. Thus, the relevant rate of 
fusion is not proportional to the square of the number of vortices but to the first power, since each 
vortex can only merge with those vortices at its perimeter that do not evolve significantly over 
long period of times. In other words, there is not an efficient mixing and hence the mean-field 
approximation underlying eq.s ( |l5|) and (|l8|) breaks down. 

As a consequence, we can write down a linear integro-differential coupled equation for the 
time evolution of the vortex populations, in close analogy to the scaling theory of linear frag- 
mentation of Ref.p^, 24 1 . A remarkable feature of this approach is that the analysis can be 



carried out for arbitrary scaling kernels. 

It is convenient to define the distribution $(S, t) of vortex energies E. For time independent 
vorticities oj, the vortex radius R is obtained from its energy from the relation E ~ u'^R'^. 
The strategy is to characterise the vortex energy population ^{E,t) by the moments m^ = 
f^ l°'<j){l)dl, where we use the scaling ansatz of the vortex energy/size distribution 

$(E,t) = i0(-), (19) 

e^ e 

and e{t) is the characteristic (time-dependent) vortex energy. The exponent —2 is required by 
(approximate) energy conservation. Note that calculating a moment corresponds to taking the 
Mellin transform of </>(/). This implies that if we know m^, we are able to retrieve (/)(/) by taking 
the inverse Mellin transform. 



Adapting from Cheng and Redner |^^ , we obtain the equations of evolution for the vortex size 
population. Assuming that the fusion rate between two vortices of similar energy E is a power 
law E^ and calling b{x) the probability that the ratio between the final vortex energy and each 
of the initial fusing vortex be x, we substitute the scaling ansatz ( |T9[ ) in the linear Smoluchowsky 
aggregation equations and obtain two separate equations with a separation constant to, one for 
the time dependence of e(t) and the other for the scaling function (p. In terms of the moments 
rua of (/), we get the following recurrence equations : 

1 — a , , 

ma+/3 = to- -nia , (20) 

ivQ — i 

and 

La= f x%{x)dx , (21) 

Jo 

From the term ^ _-^ in ([20| ) , we see that if there exist a value a* such that 

La* = l , (22) 

then all moments with a > a* will become infinite. Provided reasonable analyticity conditions 
hold, it follows that the value a* is a pole of ttIq. Taking the inverse Mellin transform of rria 
then allows us to get ^(Q and, using the existence of the pole at a*, this immediately predict 
that 

0(0 ~ /-(!+'«+"*) , (23) 

One can check that for a large variety of b{x), there are solutions of (^) with complex exponents 
a*. This signals the existence of log-periodic corrections to a pure power law distribution of 
vortex sizes. Such a scenario has been documented experimentally in the reverse situation of 
fragmentation in Ref . p5|] . This, in turn, leads to a complex separation constant w, and thus to 
a complex dynamical exponent ^ for the time evolution. 

3 Methodology of the data analysis 

Detecting log-periodic structures in noisy, intermittent data is rather subtle and apparently 
innocent data manipulations, such as using the cumulative distribution rather than the data 
itself, may completely erase such structures. Furthermore, since these structures are log-peiiodic 
it is often impossible to obtain a large number of oscillations, since numerical and experimental 
data seldom cover more than a few decades. Hence, great care must be taken not only with 
respect to what quantity to analyse in an attempt to quantify the strength and frequency of a 
log-periodic signal but also how to average different realisations. Here the question of noise also 
enters crucially, since the higher the noise level the more oscillations are necessary. The next 
section will try to shed light on these technical questions and provide guidelines for the analysis 
of the turbulent data. 

If we suspect that some distribution p (x) obeys 

p (x) ex x" (24) 

and furthermore believe a to be complex, the most direct quantification is to look at the loga- 
rithmic derivative 

^^^ = a. (25) 

alnx 

There are a number of ways to estimate the derivative from data, but they more or less 

amount to the same thing, namely locally estimating the data by some function, typically a 




50 100 150 200 250 300 350 2 2.5 3 3.5 4 4.5 5 5.5 

Figure 1: First synthetic data set and it's logarithmic derivative. 
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Figure 2: The cumulative distributions of the data set in figure |l] and it's logarithmic derivative. 





Figure 3: Lomb periodograms of the derivatives in figures || (left) and ^ (right). 
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Figure 4: Second synthetic data set and it's logarithmic derivative. 
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Figure 5: The cumulative distributions of the data set in figure Q and it's logarithmic derivative. 





Figure 6: Lomb periodograms of the derivatives in figures ^ (left) and ^ (right). 
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Figure 7: Third synthetic data set and it's logarithmic derivative. 
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Figure 8: The cumulative distributions of the data set in figure |^ and it's logarithmic derivative. 





Figure 9: Lomb periodograms of the derivatives in figures |^ (left) and |8| (right). 



nth degree polynomial and then performing the derivative. Performing the derivative directly 
on discrete data amounts to drawing a line between consecutive data points and using its slope 
as the derivative. A more sophisticated Savitsky-Golav filter [^] uses polynomials of arbitrary 
degree and preserves higher moments of the data, such as the variance, which becomes suppressed 
with a moving average (hence, a moving average is unsuitable in an analysis for log-periodicity). 
However, in the following we will use the simplest numerical procedure to calculate the derivative, 
i.e., 

dlnp{x) _ lnp{xi+i) -lnp{xi-i) 

,, \x=Xi — 1 1 V'^'-'j 

a in X m Xj+i — m Xi_i 

in order to keep the numerical manipulation transparent. 

A numerical tool that we will use extensively in this paper is the so-called Lomb periodogram 



|26|. It corresponds to a harmonic analysis using a series of local fits of a cosine (with a phase) 
with some user chosen range of frequencies. The advantage of the Lomb periodogram over a 
Fast Fourier transform is that the points does not have to be equidistantly sampled, which is 
the case here. 

3.1 The Lomb periodogram and integration 

A method often used in order to investigate fluctuations in data is to fit the average behaviour 
of the data by some function which is then subtracted. This method can be justified if you only 
want to estimate the error of that specific fit. However, you cannot use it to extract oscillations 
unless you have many and if that's the case you can do much better with the data itself. The 
problem with the procedure is that the somewhat arbitrary "truncation" of the data set in both 
ends will produce a bias due to the randomness of the phases of each realisation. If some sort 
of de-trending is to be used in order to accentuate the more subtle correlations, such as log- 
periodicity, then it always makes more sense to divide the data by the leading trend rather than 
subtracting it. However, it's a non-trivial task to get an un-biased estimate of such a leading 
trend due to the same arbitrary truncation. When this is possible, such a de-trending is a very 
powerful way of visualising log-periodic oscillations. This has been clearly illustrated in the case 



of financial data [27| 



Another troublesome feature of extracting log-periodic components from data is that they 
can easily be destroyed by "innocent" data manipulation such as calculating the cumulative 
distribution. An example where this is clearly the case is shown below. We have generated 10 
data sets with 30 points each using the equation 

y(t) = t~0-5 [1 + 0.1 -cos (27r/ In t)] with / = 1.114 (27) 

where the sampling was random in the sense that the spacing between two consecutive points 
was chosen from the interval [0 : 20] with uniform probability. 

Figures (|l|,^0) show three typical realizations of the random sampling of ( p7| ) together with 
the corresponding discrete estimations of their logarithmic derivatives. Figures (0,0J^) show the 
corresponding integrals (cumulatives) of the three samples together with the discrete estimations 
of their logarithmic derivatives. Figures (0j6|,0) compare the Lomb spectrum of the direct and 
cumulative sampling presented in the previous figures. Whereas we clearly extract the log- 
periodic component in the signal and its frequency / = 1.1 with a level of significance better 



than 0.995 (see Ref.|26] for the quantification of the significance from the peak level), the signal 
has completely disappeared in the cumulative distribution. 

This not only shows that a relatively strong (10% in amplitude) log-periodic signal will be 
destroyed by an seemingly innocent summation but also that without noise only ~ 2 oscillations 
is enough to qualify log-periodicity. Similarly to the performance of a FFT which depends on 

10 



the ratio of the measured frequency over the Nyquist frequency as well as the length of the 
time series, the performance of the Lomb periodogram depends on the number of points per 
oscillation and the number of oscillations. In fact, this example confirms very well the general 
belief that noise in data can be "integrated out" and that constructing cumulatives provides a 
very efficient low-pass (1//^) filter. However, for our present task of identifying log-periodicity, 
it also shows the danger of this procedure which integrates out both the noise and the useful 
signal. 

3.2 Additive noise 

The process used in generating figures (||,§J3) is already noisy in some sense due to the random 
samplingFI We now investigate the effect of additive noise, in addition to the irregular sampling. 
We have thus added a noise term to equation (|^) and rewritten it as 

y(t) = (l + A;(0.5-ran))t"°-^[l-h0.1-cos(27r/lnt)] , with / = 1.114 , (28) 

where k is the amplitude of the noise and ran is the intrinsic random number generator of 



Fortran. In figures |10| to |1^ we see 3 examples of such synthetic data sets and their corresponding 
Lomb periodograms. Here /c = 0.1 is used corresponding to a noise level half of the amplitude 
of the log-periodic oscillations. What we see now is that the confidence interval has come down 
considerable, but the peak is still correctly positioned for all 10 data sets of which 3 arbitrary 
sets are shown here. In fact, what we see here is not very much different from the analysis of 
diffusion limited aggregation (DLA) presented in [^, |2^. The individual data set do not carry 
a very large significance, but the consistence of the analysis for the entire ensemble of data sets 
do. 

If we increase the noise by a factor of 2, i.e., to the same amplitude as that of the log- 
periodic oscillation things become very murky, but it is still possible to identify the log-periodic 
oscillations using the results for the entire ensemble. If we, as for DLA p8| , p9[ , record the 
two highest peaks for the 10 data sets, we get the distribution shown in figure |l^. It is rather 
remarkably that, with the Lomb periodogram, we can retrieve log-periodicity with such quality 
in data which has a noise amplitude of the same size as that of the log-periodicity. In the next 
section, we shall see a more elegant way of retrieving log-periodicity in noisy data. 

3.3 The method of canonical averaging 



The method of "canonical averaging" has been introduced |3C] to reveal subtle correlations which 
would disappear in the usual ensemble averaging. We have tested a specific implementation of 
it on log-periodicity in DLA clusters and rupture time-to-failure processes |^, |2^, ^], based on 
a sample-dependent phase determined from the maximum of a sample-dependent susceptibility. 
Recall that the problem comes from the fact that the argument 27r/lnt in the cosine in 
equation (^) is defined up to a phase, since adding a phase simply corresponds to a change 
of time units. In other words, the universal physical quantity is the log- frequency / while the 



phase is expected to be sample-specific and to fluctuate from realizations to realizations |31]. If 
this problem is not addressed properly, standard averaging procedure mix measurements with 
different phases leading to a destructive interference of the log-periodic oscillations. 

We now propose a much simpler method than previously used [^l| whose principle is to "re- 
phase" the oscillations. The method is based on the fact that a Lomb periodogram is a simple 

*It is probably better to keep the vocable "noise" for processes that deteriorate the signal. In the case of 
a random sampling of n points in some interval T, frequency components which are larger than the Nyquist 
frequency of the corresponding evenly sampled points ^ can be identified. Hence, the random sampling is in 
fact improving the signal-to-noise ratio for high frequencies. 

11 
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Figure 10: First noisy data set (left) and Lomb periodograni of its logarithmic derivative 
(right). 
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Figure 11: Second noisy data set (left) and Lomb periodograni of its logarithmic derivative 
(right). 
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Figure 12: Third noisy data set(left) and Lomb periodograni of its logarithmic derivative 
(right). 
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Figure 13: Histogram of the frequency distribution for the two most significant peaks in the 
Lomb periodogram of 10 data sets with a noise amphtude of the same size as the ampHtude of 
the log-periodic oscillations. The peak in the histogram correspond to a frequency of ~ 1.12 in 
very good agreement with that of / « 1.114 used. 



cosine fit, while the phase is not playing any role. Hence, the periodogram corresponds to a 
removal of the phase from the data. This means that averaging the periodograms of different 
data sets is exactly the "re-phasing" we are looking for. In figure |l^, we show the averaged 

The two peaks 



periodogram for the 10 data sets used in the histogram shown in figure 13 



corresponds two frequencies of and 1.13 and 3.4, respectively. The first and most significant 
peak is in very good agreement with the frequency of 1.114 used in generating the data sets. 
The second peak is furthermore in excellent agreement with that of the third harmonic. We 



remind the reader that figure 14 comes from averaging periodograms of noisy data with a noise- 
amplitude of the same size as the amplitude of the log-periodic oscillations. Nevertheless, the 
log- frequency component of ~ 1.1 present in the data is unambiguously extracted as well as 
its third harmonic hardly visible in the histogram in figure 13. This illustrates not only the 



efficiency of the Lomb periodogram in retrieving periodic trends, but also that of the canonical 
averaging scheme for re-phasing the oscillations. 

This concludes the numerical tests of the tools to be used in the experimental data analysis 
that we present in the next section. 



4 Experiments on 2-d freely decaying turbulence 

The experimental data analysed here are some of those presented in |2^ (Fig.'s 6 and 7). We 
will briefiy describe the experimental techniques and the results that are used in the present 
context. 



4.1 Experimental setup 

The fiow is generated in a thin, density stratified layer of electrolyte, using an electro-magnetic 
forcing. During a short initial forcing period, an 8 x 8 array of vortices is created, with nearest 
neighbours counter-rotating. At a time defined as t = 0, the forcing is stopped, and the system 
is then decaying freely. After a vertical reorganisation of the fiow with a maximum duration of 2 
seconds, the system can be regarded as two-dimensional, except for the effect of bottom friction 
|32|, which causes an exponential decay of the total energy of the system with a well-defined 
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Figure 14: The averaged Lomb periodograms of the 10 data sets used in the histogram shown 



in figure 13. We emphasise that the amphtude of the noise was of the same size as the amphtude 
of the log-periodic oscihations. Note that the true confidence interval does not correspond to 
the values shown, since it's an average over 10 data sets. The value of the two peaks are 1.13 
and 3.4 and in excellent agreement with the 1. and 3. harmonic of the frequency / ~ 1.114 used 
in generating the data sets. 



time constant. A rescaling of time is applied under which we can consider the system as truly 
two-dimensional [p^ ]. 

The size of the experimental system is 15 cm x 15 cm, and the initial Reynolds number is 
typically 1800. The duration of the experiments is approximately 12 seconds in rescaled time 
units. 



4.2 Determination of vortex statistics 

The velocity field v{x,y) on the free surface of the fluid is determined using the Particle Image 
Velocimetry (PIV) technique described in |3^, with a resolution of 40 x 40 measurement points. 
From the velocity field, the vorticity field uj{x, y) = dxVy — dyVx is computed using a polynomial 
fit. Typically 70 velocity fields are calculated during the time of the experiment. The vorticity 
field is analysed for vortices by searching for local extrema with the absolute vorticity above a 
threshold. The results are not sensitive to the choice of threshold, and the recognition method 
has been checked with other methods. Thus we can define a number n of vortices for each 
velocity field. The position of a vortex is defined as the position of the vorticity extremum. 
This is used to extract the mean nearest-neighbour separation r of the vortices. The size of a 
vortex is defined by the connected area with an absolute vorticity above the threshold. For each 
velocity field the mean vortex radius a is computed from the mean vortex size. 



4.3 Scaling laws 

Rapidly after the forcing has been stopped, like-sign vortices start to merge. Fewer and larger 
structures are thus formed, and this continues until the energy of the system is so small that no 
further evolution of the vortices can be observed. 

We presented [^] measurements that are averages of 9 experimental realizations (in the case 
of vortex radius, only 7 realizations were used). We found that for the time period from t = 1.5 
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s to t = 10 s the quantities defined above follow the power laws 

n{t) ~ i-o.70±o.i 

a{t) ~ t0.2i±o.06 (29) 

r{t) ~ t0.38±0.08^ 

The error bars are estimated by considering the fluctuations observed between individual re- 
alizations, however without any re-phasing of the log-periodic oscillations to be extracted in 



section iA. We note that the average values for the exponents in (^) agree well with the 



scaling approach [p3 



4.4 Log-periodic analysis of the experimental data 

As explained in the previous section, three experimentally independent quantities has been 
measured as a function of time: the number of vortices, n (t), the mean vortex radius a (t) and 
the mean distance between vortices r {t). For each of these quantities, the log-periodic signatures 
present in the data will be extracted using the tools described in the section H. Specifically, 



the logarithmic derivative as defined by eq. 25 was calculated from the experimental data. 



Truncating the times series in order to avoid end-effects gave ~ 35 points for each periodogram 
to use. This corresponds well to the 30 points used in the synthetic data analysis of section ^. 

4.4.1 Number of vortices 



In figure 15, we see the number of vortices n{t) as a function of re-scaled time t (left) and its 
logarithmic derivative (right) for a specific experimental realisation. In figure |l6| (left), we see 
the corresponding periodogram of the logarithmic derivative shown in figure ^ with a peak at a 
log-frequency of ~ 4. The height of the peak is only approximately twice the height of the second 
largest which means that the confidence of the peak is not very high. Nevertheless, performing 
a simple average of the Lomb periodograms for the 9 available experimental realisations gives us 
the picture shown in figure p!^ (right). The peak has moved up slightly but is still clearly visible, 
which means that the position of a peak at a frequency of ~ 4 for the single realisation was not 
accidental, since a pure noise signal would have disappeared in an average over 9 realisations. We 
stress that a frequency of ~ 4 is approximately a factor 5 below the average sampling frequency 
of 1.6/35 ~ 22. This means that the number of points per oscillation is smaller than that of the 
synthetic data which explains the lesser performance. 

4.4.2 Mean separation of vortices 



In figure 17, we see the average separation between vortices r(t) as a function of re-scaled time t 
(left) and its logarithmic derivative (right) for an specific experimental realisation. In figure ^ 
(left), we again see the corresponding periodogram of the logarithmic derivative shown in figure 
p^ with a peak at a frequency of ~ 4. The height of the peak is again only approximately twice 
the height of the second largest. Nevertheless, performing the average of the Lomb periodograms 
for 9 experimental realisations gives us the picture shown in figure |l^ (right). The peak has 
again moved up slightly but is still clearly visible, which again means that the position of a peak 
at a frequency of ~ 4 for the single realisation was not accidental. 

4.4.3 Mean radius of vortices 



In figure 19, we see the mean radius of the vortices a{t) as a function of re-scaled time t (left) and 



its logarithmic derivative (right) for an specific experimental realisation. In figure ^ (left), we 
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Figure 15: Logarithm of the number of vortices log(n(i)) as a function of the logarithm of 
re-scaled time log (t) (left). The corresponding logarithmic derivative as a function of log (i) 
(right). 





Figure 16: Lomb periodogram of the logarithmic derivative shown in figure |T^ (left). Lomb 
periodogram of the logarithmic derivative of the number of vortices n (t) as a function of re-scaled 
time t averaged over 9 experimental realisations (right). 
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Figure 17: Logarithm of the mean separation of vortices log (r (t)) as a function of the logarithm 
of re-scaled time log (t) (left). The corresponding logarithmic derivative as a function of log (t) 
(right). 





Figure 18: Lomb periodogram of the logarithmic derivative shown in figure |T3 (left). Lomb 
periodogram of the logarithmic derivative of the mean separation of vortices r (t) as a function 
of re-scaled time t averaged over 9 experimental realisations (right). 
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Figure 19: Logarithm of the mean radius of the vortices log (a (t)) as a function of the logarithm 
of re-scaled time log (t) (left). The corresponding logarithmic derivative as a function of log (t). 





Figure 20: Lomb periodogram of the logarithmic derivative shown in figure ll^ (left). Lomb 
periodogram of the logarithmic derivative of the mean radius of the vortices a (t) as a function 
of re-scaled time t averaged over 7 experimental realisations (right). 



see the corresponding periodogram of the logarithmic derivative shown in figure ^ with a peak at 
a (log-) frequency of ~ 4. The height of the peak is again only approximately twice the height 
of the second largest. Performing the average of the Lomb periodograms for 7 experimental 



realisations available here gives us the picture shown in figure 20 (right). The peak is rather 
broad but is still clearly visible, which again signifies that the position of a peak at a frequency 
of ~ 4 for the single realisation was not accidental. 



5 Conclusion 

We have suggested that the concept of discrete time scale invariance of the underlying physics is 
more appropriate for the description of freely decaying 2-d turbulence than the usual assumption 
of a continuous time scale invariance. 

In this goal, a numerical analysis of synthetic data have demonstrated a new canonical 
averaging scheme using the Lomb periodogram as its main ingredient. We have furthermore 
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shown that seemingly innocent data manipulations completely destroy the log-periodic structures 
in these synthetic data. 

Using this new canonical averaging scheme on experimental data obtained from freely de- 
caying 2-d turbulence, we have shown that the log-periodic signatures slightly visible in the 
individual experimental realisations carry statistical significance and cannot be disregarded as 
noise when averaged over the nine available experimental realizations. It is quite reassuring 
that the analysis of the three different experimentally independent measurements all agree on a 
log-frequency of ~ 4 — 5 for the identified log-periodic oscillations. This correspond to a prefered 
scale factor A = e^'^ ~ 1.2 — 1.3 

Our study of this experiment establishes for the first time a systematic procedure for explor- 
ing the existence of log-periodic oscillations in turbulence data. In this spirit, it would be very 
interesting to analyse previous partial indications of log-periodicity in turbulent data, that can 
be found in fig. 5.1 p.58 and fig. 8.6 p.l28 of Ref. Q, fig.3.16 p. 76 of Ref.p5|, fig.lb of 



Ref.|3§ and fig. 2b of Ref. |37[. 



For a theoretical view point, we have presented a theoretical framework which makes plausible 
the presence of log-periodic oscillations in the observables. Physically, we interpret the observed 
log-periodicity as resulting from a discrete cascade of intermittent vortex coalescence processes. 
In this sense, 2-d freely decaying turbulence would not be much different from other systems in 
which log-periodicity has been documented to result from self-similar intermittent instabilities 

i- 

From the point of view of the Navier-Stokes equations, it may seem remarkable or even 
stretching credibility that log-periodicity emerges from such continuous PDE. A possible anal- 
ogy is offered by the recent discovery that the continuous nonlinear Einstein PDE's of general 
relativity in the presence of a scalar field self-interacting through gravitation may generate a 
log-periodic spectrum of black hole masses with develop according to a log-periodic self-similar 
time dynamics p8| . The mechanism might result from the existence of a limit cycle in the 
renormalisation group description of a field close to the negative density limit (in turbulence, 
could this be obtained from a negative effective viscosity?). Another possible route [^3[ is that 
scale invariant equations that present an instability at finite wave-vector k decreasing with the 
field amplitude may generate naturally a self-similar discrete spectrum of internal scales. 

We acknowledge useful comments from A. Arneodo, B. Dubrulle, U. Frisch and P. Tabeling. 
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